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1 Introduction 



With the advent of large hadron colliders, a precise knowledge of the structure of hadrons 
has become compulsory if experimental data on elementary particle physics are expected 
to be analysed and ultimately interpreted with the desired high accuracy and reliability. In 
particular, one of the goals of the LHC project is to perform precise tests of the Standard 
Model of strong, weak and electromagnetic interactions and the fundamental constituents 
of matter. Actually the LHC machine can be viewed as a gluon-gluon collider to a large 
extent and many signatures (and their backgrounds) of physics, both within and beyond 
the Standard Model, involve gluons in the initial state Therefore an accurate knowledge 
of the gluon density in protons acquires a special relevance. Currently there are three 
major groups - namely CTEQ, MRST and GRV - providing regular updates of the partonic 
structure of protons as new data and/or theoretical improvements become available ||. 

The parton distribution functions (PDFs) are phenomenologically determined by global 
analyses of a wide class of hard processes involving initial-state hadrons, making use of the 
QCD parton framework. Presently, the most precise determinations of the gluon momen- 
tum distributions in the proton basically come from data on deep-inelastic scattering (DIS), 
in particular through the analysis of the scaling violations of the structure function F 2 . 
However, this represents an indirect method since the gluon density is obtained by means 
of the QCD evolution (DGLAP) equations. On the other hand, hadron-hadron scattering 
processes with prompt photon production or jets in the final state should be (and actually 
already are) extremely adequate to probe "directly" the gluon distribution in hadrons jj. 
Such "direct" determinations can be viewed as complementary to indirect analyses, pro- 
viding independent tests of perturbative QCD with distinct systematic errors. Let us also 
remark that, in colliders like the Fermilab Tevatron or the LHC H, gluon densities are - or 
can be - probed at similar (and higher) x-values as in DIS but at significantly larger energy 
scales. 

In this paper, we examine the possibility of using heavy quarkonia inclusive hadropro- 
duction to probe the gluon density of protons, extending in a more quantitative way the 
ideas earlier presented in Refs. ||, |(| . In particular we are focusing on bottomonia production 
in proton-proton collisions at the LHC. Our claim is that, ultimately, the measurement of 
bottomonia cross sections could provide a useful consistency check of different PDF sets and 
their energy-scale evolution. 



2 Heavy quarkonia inclusive hadroproduction 

At very high transverse momentum, fragmentation mechanisms should expectedly become 
dominant in the cross section of bottomonium hadroproduction - as in any other single- 
particle production channel. Specifically, gluon fragmentation into a bb pair followed by its 
non-perturbative evolution yielding a T final state by emission of soft gluons, should play 
the leading role according to the colour-octet mechanism (COM) J7|. Such a production 
mechanism constitutes the relativistic generalization of the colour-singlet model (CSM) and 
one of the most natural explanations of the excess of heavy quarkonia production observed 
at the Tevatron P, PI pi . Moreover, the COM can be cast into a rigorous framework based 



on NRQCD ||11||, an effective theory coming from first principles. 
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Furthermore, as long as we focus on quarkonium hadroproduction at large transverse 
momentum, the factorization assumption, underlying the NRQCD description, should be 
justified by the large scale px] for bottomonia in particular the heavy quark mass is likely 
large enough to factorize the short and long distance physics in the partonic interaction 
itself, as it is well known from decay processes. Whether similar arguments can be applied to 
charmonium resonances has to be checked, for example analysing the transverse polarization 



of the resonance 12 



Therefore one can expect, on solid grounds, that one of the main production channels of 
bottomonia at the LHC should correspond to the partonic subprocess: 

9 9 -> 9* 9 (1) 
and the subsequent gluon fragmentation into a T(nS) state accompanied by light hadrons 

g*^T(nS)X s ; (n=l,2,3) (2) 

produced through the already mentioned colour-octet mechanism. 

We must clearly state that our later development relies on the hypothesis that the main 
contribution to heavy quarkonium hadroproduction at high transverse momentum comes 
from a fragmentation mechanism as in (2) whose final hadronization stage is governed by 
a single non-perturbative parameter. Strictly speaking, the colour-octet production mech- 
anism is not necessarily required; for example, another possibility is gluon fragmentation 
into a heavy resonance via a colour-singlet channel. However, this mechanism actually cor- 
responds to a higher order contribution since colour conservation and charge conjugation 
require the emission of two extra gluons in the subprocess (2) jy|. Therefore, the colour- 
octet fragmentation mechanism should play the leading role at high transverse momentum 
in bottomonium hadroproduction || if the COM turns out to be correct. In this regard, 
the different leading and subleading contributions and their relative importance will be very 
briefly reviewed in section 2.1. 

Ideally, the final-state gluon (g) in Eq.(l) will give rise to a recoiling jet (g — > jet), 
sharing, in principle, the same transverse momentum as the heavy resonance - neglecting 
initial and final state gluon radiation. On the other hand, the leptonic decay of the resonance 
is obviously the best way to have a clean signal of its formation among a huge hadronic 
background. Of particular interest is the muonic channel, since such muons would very likely 
pass the first level trigger - consisting (in ATLAS) of a muon with transverse momentum 
larger than 6 GeV and (absolute) pseudo-rapidity less than 2.5 [|IJ. 

Hence, events would topologically consist of an almost isolated muon pair from the decay 
of the heavy resonance recoiling against a single jet, therby providing a suitable tag for the 
production mechanism represented in Eqs.(l-2). Indeed this topology is the expected one 
since the mass difference between the intermediate coloured and final states should be quite 
small, allowing the emission of at most a few light hadrons via soft gluon radiation - denoted 
by X s - accompanying the resonance at the final hadronization stage. Moreover, in order 
to remove possible background, events with (more than one) mini-jets should be discarded 
from the selected sample. 
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jet (X h ) 



Figure 1: Sketch of a parton- level process leading to a resonance recoiling against a gluon yielding 
a final-state jet. Notice that if the two interacting partons had different fractions x of the total 
momenta of the colliding protons, the event would not show a back-to-back topology between the 
dimuon and the recoiling jet. Besides, initial and final state radiation can affect the event topology 
similarly. We consider an event compatible with a back-to-back topology if the rapidity difference 
in absolute value of the final-state resonance and jet is less than a predetermined value set by 



2.1 Differential cross section 

In this Section we reproduce the main points presented in Ref. for the sake of clarity in 
our later more technical proposal and critical discussion. 

According to the collinear approximation, employed in our analysis, the triple differential 
cross section for the inclusive production process pp —>■ TX can be written as 

- — j ° = 2p T x a x b f a/p (x a , Q 2 )f b /p{x b , Q 2 ) d(Jah ^ ) (3) 
dy r dy jet dp T dt 

where px denotes the transverse momentum of both the resonance and the jet; yr and yj et 
represent the corresponding rapidities in the Lab frame, respectively; f a ,b/p(x a ,b, Q 2 ) denotes 
the parton-a, b density in the proton at Feynman x a , b and typical hard interaction scale Q 2 . 
(Let us note that all factorization, fragmentation and renormalization scales have been taken 
equal for simplicity.) 

The partonic differential cross section can be written as: 

consisting of a short- distance (and calculable) part, i.e. A(ab— >[bb] n X h ) where X h denotes 
the hadron bulk from the hard interaction, as depicted in Fig.l, and a long-distance piece 
(O^) which can be identified as a NRQCD matrix element The barred summation 

represents an average over initial spin and colour for the short-distance partonic process, 
and there is another summation running over channels corresponding to different quantum 
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numbers of the intermediate state (labelled by n) evolving into final heavy quarkonium. 
Using the spectroscopic notation for the [bb] n state: 2S+1 L^j' 8 \ where the superscript (1,8) 
stands for either a colour- singlet or a colour-octet state, the different contributions are: 
^Sq , 3 S[ 8 \ ^Sq^, 3 Pj 8 \ ... The relative importance of each contribution is basically given by 
its pr dependence, the order of the strong coupling constant a s and the so-called velocity 
scaling rules of NRQCD JT]] (a nice overview can be found in [TJ[]). In this work, we shall 



concentrate only on direct production at high-pr, (dominated by the 3 Si contribution) 
without feedown from Xbj states. 

The basic kinematics of the partonic subprocess (see Figure 1) is determined by the 
transverse momentum p? and the rapidities of the resonance yr and of the recoiling jet yj et , 
respectively []; therefore yy + yj et is twice the rapidity of the partonic centre-of-mass system 
yo (i.e. yo = (yy + yj et )/2) and 5y = yr — yjet denotes the rapidity difference. 

At large transverse momentum one can write: 

x a , b ~ 2PT e — cosh^ (5) 

where the measured transverse momentum of the resonance pr is assumed to coincide at 
leading order with the transverse energy of the jet. Eq.(5) can be seen as a particular ap- 
plication of well-known formulae for dijet hadroproduction (see ]Tj| and references therein). 

At relatively high p? the dominant partonic subprocess for T hadroproduction should 
be the gluon-gluon interaction |J. Thereby, the following combination appears in the cross 
section of Eq.(3) particularizing to the case of gluon colliding partons a and b: f g / p (x,Q 2 ) 
= g(x,Q 2 ), 

x a x b g(x a ,Q 2 ) g(x b ,Q 2 ) 

calculated at the common scale Q 2 , basically set by p\ as we shall see. 

In particular, requiring y = (i.e. \yx\ = \yj e t\ but with opposite signs) implies consid- 
ering only those values x = x a = x b (within the uncertainty interval to be discussed below) 
given by 

x = xt coshy T (6) 



with xt = 2pT/ys. Thus the theoretical and experimental analysis simplifies notably: 
only the transverse momentum and rapidity of the T resonance enter into Eq.(6), both 
determined from the experimental measurement of the transverse momenta and pseudo- 
rapidities of the muons from its decay. 

On the other hand, the Feynman x can be known with a precision limited by the rapidity 
uncertainty of the resonance; indeed, for a fixed px value, the typical uncertainty Ax is given 
by! 

< Vent (7) 

X 



1 According to energy /momentum leading-order balance of the parton-level interaction, we assume that the rapid- 
ity of the parent gluon equals the jet rapidity which can be obtained from its measured pseudo-rapidity r\j at — yjet 
(for technical details about the measurement of the jet pseudorapidity in the hadronic calorimeter we refer the reader 
to |). 
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where y cu t represents the half-width of the allowed rapidity interval for the resonance, cor- 
responding to a back-to-back topology as shown in Fig.l 0. 
Then we can write as a first approximation, 

d3(X - 2p T x 2 g(x, Q 2 ) 2 da ^ Q2) ( 8 ) 



dyrdyjetdpT ' '' dt 

On the other hand we can choose the typical hard scale as the momentum transfer 
Q 2 = —t, where t stands for the Mandelstam variable of the partonic subprocess. This 
quantity can be estimated by means of the expression 

Q 2 = 2p\ cosh 2 ^ ^1 — tanh^^ = 2p 2 T cosh 2 y r (1 — tanhy T ) (9) 

the last step coming since 5y = 2y r as we have assumed y = 0. 

Let us remark that x and Q 2 are not completely independent since the centre-of-mass 
energy of the partonic interaction s = x 2 s can be expressed as s = Ap\. cosh 2 yx, neglecting 
masses |. Therefore the following relation is satisfied 

Q 2 = |x 2 (1-tanhyx) (10) 

Moreover, keeping x fixed and averaging over the full rapidity range, we get 

(Q 2 ) = \ * 2 (ii) 

In fact, in our later development we will integrate cross sections over the full available 
rapidity range and this average value for Q 2 will be extensively used. 



3 Probing the gluon density in the proton 

In Ref. 0, we proposed to test different parameterizations of the gluon distribution g(x, Q 2 ) 
in the proton by studying the ratios 

x 2 g(xj,Q 2 ) 2 = / dagg/dij \ ^ /p Ti \ / d 3 a/ dy r dy jet dp T j \ ^ 
x 2 g(xi,Q 2 ) 2 Kddgg/dijJ \p T jJ \d 3 a/dy T dy je tdp T iJ Vjet=-yr 

corresponding to several of the parameterization set under consideration. Subin- 

dices i,j denote discrete values or points of the x variable associated to bins whose relative 
widths are set by Eq.(7); Q 2 stands for the (average) transverse scale according to Eq.(ll) for 

2 We have assumed that the uncertainty in the determination of the Feynman x value provided by Eq.(7) deter- 
mines the binning of this variable. Thus, the x-bin size is ultimately set by the value of ycut, used to impose the 
back-to-back condition on the events (i.e. the measured rapidities of the resonance and jet should differ in absolute 
value by less that twice y C ut)- We have explicitly checked that the difference, event by event, between the "observed" 
x a .b (through Eq.(5)) after applying ATLFAST jet and T reconstruction [jjj], and the corresponding values selected 
by PYTHIA in the hard partonic interaction, is of the order of 10%, compatible with our choice y cut = 0.2. (we are 
indebted to F. Camarena for his valuable help in this calculation). 

3 We have neglected masses only in the estimate of Q 2 , not in the generation of events with PYTHIA. Actually 
for py^A/- 2 -, this leads to no significant difference in the former case. 
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x = X{. Actually the ratio given by Eq.(12) involves the so-called effective gluon distributions 



nj, since quark distributions can contribute to some extent via the partonic scattering 
channel gq—>-g*q, and to a much lower extent via the channel qq-^g*g @. 

It is important to note that the number of independent x values is basically limited by 
Ax (i.e. by y cu t), and hence the number of available Xi,Xj pairs in the ratio (12); therefore 
this condition cannot be released too much. At the same time, the foreseen statistics does 
not allow one to decrease y cut too much either, and a compromise should be reached (our 
particular choice has been y cut = 0.2). 

In sum, the keypoint of our proposal is to consider the l.h.s. of Eq.(12) as an input 
corresponding to different sets of the gluon distribution for the proton, whose x dependence 
is hence assumed to be "known", and in fact would be tested. On the other hand the r.h.s. 
corresponds to an input from experimental data and some theoretical factors likely under 
control. 

For the sake of clarity let us rewrite expression (12) as: 



where 



Rtheo x Rexp (13) 



Rtheoivn , vtj , Pn , ptj ; Ql ,Q 2 j) = d jr^nf ( 1 1 > 



d&gg/dt 



J 



In the high p? limit and zero rapidity (yn = yrj — 0), 

<4(Qi) PTi 



R 



theo 



m) P 4 Ti 



it 

\2\ 



explicitly showing that a s (Q ) is entangled in the gluon density determination. Indeed, the 
determination of a s from the analysis of inclusive jet cross sections at hadron colliders can 
be viewed as deeply related - in a complementary way - to the study of the sensitivity of 
the hadronic production cross sections to different PDF sets fT5|| . 

Since at high p T only the 3 S^ channel will be considered (see next Section), the depen- 
dence on the NRQCD matrix elements does cancel in Rtheo, but there is a dependence on the 
scales Q\ and Q 2 a which should match the corresponding dependence in the left hand side. 
Let us stress that to achieve the cancellation of NRQCD matrix elements in the quotient, 
only a single channel should contribute dominantly to the cross section. This fact exludes 
including subleading contributions, preventing us, in principle, from considering any pos- 
sible check of next-to-leading PDFs according to our method, and restricting the test to 
leading-order gluon distributions. 

On the other hand the experimental input reads as the ratio 

D / x (PTi\ fd 3 a/dyrdy jet dpTj\ 

Rexp{PTi,PTj; Vr)= X 

\PTjJ \d i a/dy T dy jet dp Ti J 
which can be obtained directly from experimental data f\ 

4 Let us note that in the condition yj et = —yr imposed to the events - within the uncertainty interval - could 
be removed, that is not requiring a back-to-back topology. This would amount to extend the study of the ratio (12) 
over the region of the plane x a ,Xb beyond the diagonal region x a — Xb, thereby allowing additional cross-checks of 
the PDF under scrutiny; however, we have not included this study in the present paper. 



(15) 

Vjet=-yr 
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Figure 2: Predicted T(15) + T(2S) + T(3S) weighted contributions to bottomonia inclusive pro- 
duction at the LHC in the rapidity interval \yr\ < 0.2 and px > 10 GeV. Only direct production 
is considered (no feed-down from Xbj states). Left panel: differential cross section; right panel: 
integrated cross section (see Ref. ||). Let us note that once the colour-octet MEs have been 
determined from the Tevatron data, using either CTEQ4L or CTEQ5L as the input PDF, the 
respective predictions for the LHC are very similar [|T^|, 



3.1 Efficiencies, statistics and expected accuracy 

From an experimental point of view, one of the advantages of heavy quarkonium physics 
(bottomonium in particular) at the LHC is its clean (and self-triggering) signal through 
the muonic decay channel: T(nS) — > fi + n~ allowing a full kinematic reconstruction. How- 
ever, it may happen that the discrimination among the different T(nS) states via mass 
reconstruction could become a difficult task at very high px because of the uncertainty on 
the measurement of quite large momenta of muons. Indeed, above px — 50 GeV for each 
individual muon momenta, the uncertainty associated to the reconstruted mass from two 
muons makes impossible the discrimination below the 1 GeV level [I]. Nevertheless, since 
we are proposing to study ratios of cross sections, we can consider the overall T(nS) direct 
production, without separating the different bottomonia sources - the weighted matrix ele- 
ment cancelling in the quotient (14) if the mass difference between the different T states is 
neglected [§]. 

In order to assess the foreseen statistics and efficiency factors, we have employed the 
PYTHIA f\ event generator ||16[ with the colour-octet production mechanism implemented 



as a new code in the Monte Carlo programme (a detailed account of the implementation, 
results and discussion can be found in Ref. B). For the purpose of illustration we plot 
in Fig.2 the combined production rate at px > 10 GeV and \yr\ < 0.2 for the upper and 

5 The quarkonium (charmonium) sector regarding 2^2 production channels has remained unchanged in PYTHIA 
from the older version 5.7 - used in our earlier studies - to the current one 6.2 - leading to the same results in all 
cases for high pr- 
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lower values of the colour-octet matrix element obtained in |J from previous fits to CDF 
experimental data of Run IB at the Tevatron. 

Note that despite the respective muonic branching ratios are already included in the plot, 
actual measurements require taking into account efficiencies for triggering, reconstruction 
and identification of particles and jets from detectors. In our study we have not performed 
a full simulation of the detector effects, although we have benefitted from different studies 
focusing on leptons and jets separately 0. Since we are addressing the feasibility of this 
proposal, rather than searching for precise predictions, this approximation should suffice. 

Firstly, the reconstruction of the T(nS) mass from the pair implies a reduction 

factor of about 80% by applying the same cuts as for the J/ip |IJ], i.e. a [—3a, +3cx] symmetric 
window around the nominal mass. These figures do not include however the muon trigger 
and identification efficiencies which altogether roughly amount to 85% for the triggering 
muon and 95% for its partner, respectively. Therefore, the overall reconstruction factor can 
be estimated as 

e M+M _ = 0.8 x 0.85 x 0.95 ~ 0.65 

Moreover, jet reconstruction efficiency, e jet , amounts on the average to about 75% 
in the px range under consideration, using an algorithm with a jet cone radius AR = 
^An 2 + A<f)' 2 = 0.7 fl8" l- On the other hand, there is an additional reduction factor, de- 



noted as e y , due to the requirement of a back-to-back topology compatible with the con- 
straint given by Eq.(7), which sets the expected accuracy in the determination of x\ choosing 
Ucut = 0.2 as already mentioned, we found (see Ref. || for more details) that e y ~ 35%. 
In sum, the combined efficiency and reduction factor to be applied is 

(-tot = V>~ x (jet x ( y = 0.65 x 0.75 x 0.35 ~ 0.17 

Assuming an integrated luminosity of 10 fb _1 , corresponding to one year running (10 7 s) 
of LHC at "low" luminosity (10 33 cm~ 2 s _1 ) and sweeping all the rapidity range [-2.5,2.5] by 
steps of 0.4 rapidity units - according to our choice y cut = 0.2 - we find from our simulation 
an expected number of about 10,000 events, corresponding to the highest Feynman x bin 
around the value of 0.014 considered in our analysis. Multiplying this number by the total 
efficiency factor e tot , we get the estimate: 

~ 1, 700 events per year run 

Hence after three years of LHC running at low luminosity f], the accumulated statistics at 
the largest x bin would amount to about 5, 000 events, basically fixing the typical accuracy 
of the order of ~ 2.5%, which essentially sets the "discrimination level" to be applied to 
different sets of PDFs, as we shall see in the next section. 



4 Testing the x-shape of different PDFs 

As already pointed out, the l.h.s. of Eq.(12) can be considered as an input from a given 
parameterization set of the gluon distribution of the proton; in Fig. 3.b we show the re- 
shape for different choices (CTEQ, MRST, GRV) |B|. All distributions were normalized 

6 We are considering the possibility of extending our proposal to the high luminosity regime at the LHC coping 
with pile-up effects, as well as using the electronic decay channel of the T resonances, too. This would mean to 
increase considerably the foreseen statistics allowing to raise the x range available in our study. 
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Figure 3: a) Xig(xi,Q 2 ) values for different leading log (LL) PDFs [19|. b) Ratios 
x 2 g 2 (xi,Q 2 )/x 2 g 2 (x^,Q 2 ), i.e. values of x 2 g 2 {xi, Qf) normalized to the rightmost point (i = 5); Q 2 
represents here an average value over the whole rapidity according to Eq.(ll) but dependent on 



to the rightmost point of Fig. 3. a, i.e. corresponding to x§ = 0.014. (Although arbitrary in 
principle, this particular choice seems quite appropriate since all PDFs become close each 
other near this point.) The points shown in the plot correspond to different values of the 
x variable whose bin widths follow the requirement of Eq. (7) with y cut = 0.2; Q 2 was 
selected as the average value of Eq. (11) for each x value. For the sake of clarity, let us 
stress that the plots in Fig. 3 have a "general validity", i.e. can be obtained from a any 
PDF parameterization set under the above-mentioned conditions, without any reference to a 
specific experiment; the applicability to a particular case depends crucially on the expected 
number of events. 

In this regard we show in Fig. 4 the maximum difference (in %) obtained from the 
curve of Fig.3.b, amounting to about 35% at most over the whole x interval. Therefore 
the foreseen precision level of 2.5% (obtained in the previous Section) likely should permit 
discriminating easily among those LL PDFs at the LHC Finally, let us note that, according 
to our remark in Section 2, we have not included next-to-leading order PDFs in our analysis 
since this would require, for consistency, further subleading contributions even at high p?, 
and therefore new NRQCD matrix elements not completely cancelling in the ratio (14). 

On the other hand, it might be possible to perform an overall consistency check on bot- 
tomonium production by considering the integrated cross section above a lower pTmin cutoff, 
say 20 GeV, for different parton distributions. This is in accordance with the suggestion 
made in of looking at the uncertainties on physical quantities rather than obtaining 

errors associated to the PDFs. In particular, expected differences of the theoretical calcula- 
tions on the integrated cross sections above PTmin = 20 GeV for bottomonia hadropro duct ion 
at the LHC, using different parameterization sets, is of order 40% respectively []. 

7 For a lower pr-min the difference becomes even larger; for instance, it could reach a factor 2 for pr > 1 GeV! 
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Figure 4: Differences (in %) between current "normalized" LL PDFs corresponding to the upper 
(CTEQ) and lower (MRST) values of Fig. 3.b. A foreseen precision of ~ 2.5% should suffice to 
discriminate among them over the whole x range shown in the plot. 

5 Conclusions and last remarks 

With the advent of high energy, high luminosity proton colliders, it is universally recognized 
the importance of knowing the uncertainties of the parton distributions to perform accurate 
measurements of the Standard Model parameters and carry out a tantalizing discovery 
programme on new physics. Keeping this interest in mind, "direct" determinations of the 
gluon density by means of distinct processes should permit the computation of systematic 
and theoretical errors, as well as putting stringent constraints to global fits where DIS data, 
however, are still playing a central role. 

To this aim and in view of the foreseen rates of bottomonia production at relatively large 
transverse momentum in the LHC 0, we propose testing the a;-shape of the gluon density 
in protons. The proposal could be similarly applied to forthcoming experimental data from 
the high-luminosity run II at the Tevatron. 

If gluon fragmentation is confirmed as the dominant production bottomonia mechanism 
at high transverse momentum, the topology of events would consist of an almost isolated 
lepton pair coming from the resonance decay, recoiling against a jet. This represents basically 
two advantages: 

• There is a well reconstructed direction defined by the lepton pair which can be assigned 
to the fragmenting gluon in the final state of the hard interaction, thereby helping in 
the search for the partner jet with a back-to-back topology, tagging altogether the 
event . 

However at small/moderate pr the colour-singlet mechanism would be largely dominant with quite more parameters 
entering the computation - such as the wave functions (and derivatives) of bottomonia states, cascade branching 
ratios M - which would make more complex the physical interpretation of any discrepancy with the experimental 
result. 
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• According to the colour-octet production mechanism, the lepton pair should provide, 
by means of its measured rapidity and transverse momentum, a clean kinematic infor- 
mation about the partonic subprocess itself (i.e. the x- values of the colliding partons). 

From inspection of the Figs. (3-4) and the expected discrimination level (~ 2.5%) cor- 
responding to the foreseen statistics collected after three years of data-taking at the LHC 
(at "low" luminosity), we concluded that a clear discrimination between different sets of LL 
PDFs should be feasible. In this regard, let us stress that Monte Carlo simulations in high 
energy physics have definitely become an ordinary - as well as necessary - tool nowadays for 
the experimental analysis of data, in particular to disentangle interesting signals from back- 
ground. General-purpose event generators implementing LL PDFs as default options (like 
PYTHIA) are widely used by the scientific community involved in the front-line analysis 
of experimental data. For instance, cross section measurements often require Monte Carlo 
estimates of the fraction of unobserved events due to detector acceptance and efficiency 
limits; an equivalent argument applies to background estimates. A knowledge, as precise as 
possible, of the dependence of the theoretical predictions on the choice of a particular PDF 
is mandatory. 

Then, we have concluded that our proposal could be useful as a consistency check of 
different gluon density parameterizations for the proton, of interest for other physical pro- 
cesses involving gluons in the initial state, too. In particular, exciting searches for extra 
dimensions have being recently carried out at the Tevatron [^TJ and are foreseen at the 



LHC [^] where effective Planck scales |23j of the order of the TeV can be tested. Relevant 
processes in the latter case are gg — > Gkk9, 91 — > GkkQ, ■■■ where Gkk denotes a Kaluza- 
Klein graviton, leading to a final-state single jet and missing energy. The channel given 
by Eqs.(l-2) should also be helpful in the study of the detector response to a monojet and 
missing transverse energy, taking into account in the latter case the information provided 
by the recoiling muonic pair. 

On the other hand, parton distributions unintegrated over the parton transverse momen- 
tum are increasingly becoming employed in the analysis of physical processes initiated by 
hadrons. Indeed, one of their advantages is that they correspond to the quantities entering 
Feynman diagrams allowing for true kinematics even at leading order. In our particular case, 
this would demand to modify the cross section given in Eq.(3). However, the theoretical 
situation regarding heavy quarkonium production is still controversial since there are impor- 
tant discrepancies between the collinear approximation and the so-called factorization 

m. 



A final remark is in order. Although the universality of the colour-octet matrix elements 



is not definitely well established |13| , one could expect that the matrix elements obtained in 
essentially the same kind of hadronic processes (i.e. hadroproduction at the Tevatron) should 
become reliable enough once used at LHC energies, under the same theoretical inputs (e.g. 
the bottom mass, choice for the factorization and fragmentation scales, etc). Therefore, 
if a precise normalization of the bottomonia cross section is provided (using forthcoming 
data from high-luminosity Run II at the Tevatron) the possibility exists of probing different 
PDF sets by means of an unfolding procedure from heavy quarkonia hadroproduction at 
the LHC. 
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